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Abstract 

The transition matrix elements for the 0+ 0^ double beta decays are calculated for 
^^Ca, ™Ge, ^"^Se, ^°°Mo, ^^^Te and ^^°Te nuclei, using a ^-interaction. As a guide, to fix 
the particle-particle interaction strengths, we exploit the fact that the missing symmetries 
of the mean field approximation are restored in the random phase approximation by the 
residual interaction. Thus, the T = 1, S" = and T = 0, S* = 1 coupling strengths 
have been estimated by invoking the partial restoration of the isospin and Wigner SU(4) 
symmetries, respectively. When this recipe is strictly applied, the calculation is consistent 
with the experimental limit for the 2v lifetime of ^^Ca and it also correctly reproduces 
the 2v lifetime of ^^Se. In this way, however, the two-neutrino matrix elements for the 
remaining nuclei are either underestimated (for ™Ge and ^^'^Mo) or overestimated (for ^"^^Te 
and ^^°Te) approximately by a factor of 3. With a comparatively small variation (< 10%) 
of the spin-triplet parameter, near the value suggested by the SU(4) symmetry, it is possible 
to reproduce the measured TIJ2 in all the cases. The upper limit for the effective neutrino 
mass, as obtained from the theoretical estimates of Oz/ matrix elements, is < >= 1 eV. 
The dependence of the nuclear matrix elements on the size of the configuration space has 
been also analyzed. 
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1 Introduction 



The odd-odd isobar, within the isobaric triplet {N,Z), {N - l,Z + 1), {N -2,Z + 2), 
generally has a higher mass than its even-even neighbors due to the pairing energy. As such, 
two consecutive (3 decays are energetically forbidden. Yet, the initial nucleus (A^, Z) may 
decay to the final nucleus (A^ — 2, Z + 2) through virtual excitations of the intermediate 
nucleus {N — 1, Z + 1), causing a double beta {[3(3) decay process. In case the lepton number 
is strictly conserved the neutrino is a Dirac fermion [v ^ v) and the two-neutrino mode 
{Pi32u) is the only possible mode of disintegration. On the other hand, if this conservation is 
violated, the neutrino is a Majorana particle {u = u) and the neutrinoless double beta decay 
iPPou) also can occur. 

Both the Pl32u and /3/3ov-decay processes have attracted much attention during the last 
decade. This is mainly because the neutrinoless decay mode constitutes the most critical 
touchstone for various gauge models that go beyond the standard SU{2) l x U{1) gauge model 
of the electroweak interaction. The (3(3^^ decay rate depends on several unknown parameters 
(neutrino mass, majoron coupling, the coupling constants of the right-handed components of 
the weak Hamiltonian, etc.) and the only way to put these in evidence is by having sufficient 
command over the nuclear structure. It is precisely at this point that the I3(32y decay mode 
is important. A comparison between experiment and theory for the /5/?2!/-decay, provides a 
measure of confidence one may have in the nuclear wave functions employed for extracting 
the unknown parameters from (3(3Qy lifetime measurements. 

In recent years, following the work of Vogel and Zirnbauer [||], the proton- neutron (pn) 
quasiparticle random phase approximation (QRPA) has turned out to be the most popular 
model for calculating the nuclear wave functions involved in the /5/5-transitions [0, ^, ^ |^, 
I, 0, I, I, |10|, 0, ig, |13|, 0, |15|, |16|, 0, m, ig, P, |2[. The QRPA calculations have 



shown that the ground state correlations (GSC), induced by the residual pn interaction 
in the particle-particle (PP) channel, play an essential role in quenching the I3l32v decay 
probabilities. Regarding the sensitivity of /?/3ov nuclear moments, M-ou, to the PP interaction, 
no consensus has been reached so far. The results obtained by the Pasadena group , with 
a residual 5 interaction, seem to show that M.qi, are rather sensitive to this force and the 
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^^Ge lifetime measurement (^1/2 > 1-4 x lO^'^yr) yields the upper bound on the effective 
neutrino mass, < rrii, >, varying between 4.4 and 10 eV (see table |1|). On the other hand, 
the calculations done by Tiibingen and Heidelberg groups with a realistic interaction (the 
G-matrix of the Paris and Bonn potentials) [|, H, |18|, suggest that the A^oi/ moments 



are to some extent insensitive to the PP channel, which results in a smaller neutrino mass 



with more stringent upper limit (between 2.3 eV and 3.1 eV) from the same datum |2^. The 
difference between the three sets of calculations lies in the importance attached to the role 
of destructive interference between the allowed (L = 0) and the forbidden (L > 0) virtual 
transitions. There is a general agreement that, for the physically acceptable values of PP 
coupling strength, L = contributions are relatively small and those from L > processes 
are significant. But, whereas the Pasadena group believes that the allowed matrix elements 
are strongly canceled by the forbidden ones, the other two groups claim the cancellation 
effect to be of minor importance. The exceptions are the results for the Aiou moments in 
^^Ca, ^°°Mo and ^^^Te: for the first two nuclei very small nuclear matrix elements have been 
obtained by Heidelberg group, suggesting that the L = matrix elements are sizable, while 
for the last nucleus all three groups obtain similar results. 

An important question in the QRPA calculations is, how to fix the PP interaction 
strengths within the S = 0, T = 1 and S = 1, T = channels? Several attempts have 
been made to calibrate the last one using the experimental data for individual Gamow- Teller 
(GT) positron decays @, 0|. The weak point of this procedure is that the distribution 
of the (3~^ strength among low-lying states in odd-odd nuclei is certainly affected by the 
charge-conserving vibrations, which are not included in the QRPA. Q In the present work 
we resort to the restoration of the isospin and spin-isospin SU(4) symmetries to gauge, re- 
spectively, the 5* = 0, T = 1 and 5* = 1, T = strengths in the PP channel 0. Unlike 
the method mentioned above, this procedure involves the total Fermi (F) and GT strengths, 
which dependent of the charge-conserving vibrations only very weakly. We are aware, how- 



^The single beta transitions ^'^'^Tc ^'^'^Mo and ^'^°Tc have been discussed recently in the 

standard QRPA However, from the analysis of the structure of the triplet of low- lying states in ^°^Afo 
(cf. ref. |3^), it can be inferred that the collective degrees of freedom should play a very important role in 
these decays. 
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ever, that the SU(4) symmetry is badly broken in heavier nuclei like those considered here. 
As such, before proceeding further, it is necessary to specify what we mean by reconstruction 
of isospin and SU(4) symmetries in the context of our calculation. 

For a system with N ^ the isospin and spin-isospin symmetries are violated in the 
mean field approximation, even if the nuclear hamiltonian commutes with the corresponding 
excitation operators (t=p and o"t=p). But, Thouless [^] and Brown have shown that 
when a non-dynamical violation occurs in the Hartree-Fock (HF) solution, the RPA induced 
GSC can be invoked to restore the symmetry |3^. There are subtleties involved in 

the restoration mechanism: the GSC are not put in evidence explicitly, but only implicitly 
via their effects on the one-body moments (S^^ between the ground state and the excited 
states. Besides, for the F excitations and when the isospin non-conserving forces are absent, 
a self-consistent inclusion of the GSC leads to the following: 

i) all the (3~ strength is concentrated in the isobaric analog state (IAS), and 

ii) the j3^ spectrum, which in RPA can be viewed as an extension of the /3~ spectrum to 
negative energies, is totally quenched. 

The extent to which the above conditions are fulfilled may be taken as a measure of the 
isospin symmetry restoration. Lee has shown that even when the isospin symmetry is 
dynamically broken by the Coulomb force, the RPA induced GSC result in a near restoration 
of the isospin symmetry. Besides being spontaneously broken by the HF approximation, the 
SU(4) symmetry is also dynamically torn down by the spin-orbit field and the supermultiplet- 
destroying residual interactions. But, as argued below, the last two effects have a tendency 
to cancel each other. 

In the mean field approximation (and due to the spin-orbit splitting) the (Tt_ operator 
has roughly equal strengths for the spin flip as compared to non-spin flip transitions. Charge 
exchange reactions |^6|, ^ ^ have revealed, still, that the spin-isospin residual forces can to 



some extent overcome this SU(4) symmetry breaking by transferring a substantial amount 

of the non-spin flip strength to higher energies (and build up in this way the GT resonance). 

^We may point out that it is not possible to obtain similar results in the Tamm-Dancoff approximation 
(TDA) where the GSC are neglected. As a matter of fact within the TDA the /3~ strength is always 
fragmented and the perturbed strength remains equal to its unperturbed value. 
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In fact, the observed mean energy differences between the GT and F resonances can be 



accounted for by the relation |39 



GT 



2QA-^/^ - 18.5 



N -Z 
A 



MeV, 



and the displacement of the GT resonance towards the IAS with increasing N — Z may be 
interpreted as the effect of the residual interaction |T|. As a matter of fact, the first 

term has the mass dependence of the mean spin-orbit splitting, A/^, and within a schematic 
TDA for the S interaction (that we use in the present work) one gets 



GT 



Ef 



A 



Is 



N - Z 
2A 



MeV, 



where Vs and Vt are, respectively, the singlet and the triplet coupling constants. It should be 
noticed that the GSC are likely to alter this result very little. But, within the RPA the crt+ 
transition strength is strongly quenched and the GT resonance is somewhat narrowed, as 
compared with the TDA results [Q. As such the global effect of the pn residual interaction 
on the GT strength is qualitatively similar to the corresponding effect on the F strength, 
and we say that the SU(4) symmetry is partially restored. Finally a word of caution, the 
closeness of GT resonance to IAS does not guarantee a complete reconstruction of the SU(4) 
symmetry. For example, the observed GT strength in ^osp^ jg j^-^ ^ single continuous 
peak located at the energy of the IAS, but simultaneously the GT resonance is rather broad 
(fti4Me1/; see ref. 

In the present work we perform a detailed analysis of the /3/3oiy decay rates for several 

nuclei, always focussing our attention on the < irii, > term (mass mechanism). The term 

containing the right-handed leptonic currents (RHC mechanism) is not considered. Besides 

that, no reference is made to the neutrinoless decays with single PPq^^b, and 

■^The problem of reconstruction of the Wigner SU(4) symmetry and its relationship with the proton- 
neutron interaction, has been discussed earlier in the context of the GT transitions by Fujita and Ikeda 
and by Bohr and Mottelson [Q. Recently Bernabeu et al. |Q have shown that treating the SU(4) 
symmetry-breaking part of the interaction in perturbation theory leads to M2v moments of correct order of 
magnitude. Several other manifestations of the SU(4) symmetry in nuclei have also been reported |45[ 
1811910 
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double PPou,2B, majoron emissions. The PPou,b decay has been ruled out by the mea- 



surement of the Z width |^5[, as shown by Gonzalez-Garcia and Nir ||56|, whereas no data is 
available for the PPoij,2B- We may note that the effective couplings between neutrinos and 
the majorons are easily estimated once the nuclear matrix elements are known. Presently 
our main aim is to set up an upper limit on the neutrino mass from the Oz/-decay analysis. 
Since a good understanding of the Oz/-decay requires that we understand the 2i/-decay really 
well, the PP2U decay rates have been reviewed as well. 

The paper is organized as follows. In section |^ we sketch the main formalism and list 
various formulae needed to compute the half lives for the 0^ O"*" transitions. Section ^ deals 
with the discussion of some important features of nuclear matrix elements, the comparison 
of numerical results with the experimental data, and the present limits on the neutrino mass. 
Finally, the conclusions are drawn in section |^. The evaluation of radial integrals, involved 
in the Oz/-decay, is shortly reviewed in appendix A. 

2 Formalism 



To apply the Horie-Sasaki method l^^l' write down the Fourier-Bessel expansion of Ou 
transition amplitude as 

-^0. = ^ E S [ dqYiq;uj^j.){0^\e'^'^a^t+\arM) ® {a.r M\e-'^''^alt+\Ot) , 

with 

,,21 



7rq{q + uJaj^) 
being the neutrino potential, and 

i^aj^ = EaJ^ — [Ei + Ef)/2. 

Here R is the nuclear radius introduced to make Aiou dimensionless, the symbol stands 
for the scalar product of spin matrices as defined in ref. ES] (with 0"^ = 1 for S = 0, = a 



for S = 1 and S = y/2S + 1), the summation goes over all virtual states | aJ'^M), and Ei 
Eaj^ and Ej are the energies of the initial, intermediate and final state, respectively. 
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After integrating over angles and summing up over M, we get 



X (Ot I |C+(gri; LSJ) \ |a J") (a J"| |C+(gr2; LSJ) 1 10+), 



where the quantities 

0+(gr; LSJ) = jL{qr)iYL ® a')j t. 

are one-body operators. 



Within the QRPA formulation presented in ref. |]IU[ the energies cOaj^ are solutions of 



the QRPA equation, and the transition matrix elements are given by 

{ar\\O+{qr;LSJ)\\0t) = -J2{p\\0{qr; LSJ)\\n)A+{pn; aJ^ 

pn 

{0^\\O4qr;LSJ)\\ar) = -J2{p\\0{qr;LSJ)\\n)A*_{pn;ar) 

pn 

where the reduced pn form factors are 

V^(p||0(gr;L5J)||n) = ^^^^/^.(g), 

with the radial part 

roo 

Rpn{(l) = Un{r)up{r)jL{qr)r'^dr, 
u{t) being the single-particle radial wave functions, and the angular part 



-P 2 JP 



The amplitudes A-|-(pn; aJ'^) are defined as: 



A+(pn; aJ) = ^PpPn [UpVnXpn-^aJ + VpUnYp^.^j] , 

A_(pn; a J) = [vpU^Xp^-aj + UpV^Yp^.^ij] , 
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where the unbarred (barred) quantities indicate that the quasiparticles are defined with 
respect to the initial (final) nucleus; = + v^, = + v^, and all the remaining 
notation has the standard meaning |TT| , |22| . As such, we obtain 

moAr) = E A+(pn;aJ")A*_(p'n';aJ^) 



with 



LS 



POO 

7^^(pnp'n';u;„J.) = ^ dqq\{q;uJc,j^)RpAq)Rp,^^,{q). 



(2) 



(3) 



The expression (|^), for the evaluation of the Ou moments within the QRPA, is much simpler, 
regarding the angular momentum recoupling, than the ones currently used in the literature 



181 , The radial integrals (y) are of the same type as the ones that appear in the 



work of Horie and Sasaki so we use their method to evaluate these, without getting 
involved with Moshinsky brackets. To show the simplicity of this procedure, a few results 
relevant to this work are reviewed in appendix A. 

In the present work we consider both the GT and Fermi 2z/ transition matrix elements 
with the corresponding amplitudes written in the form: 



M2u = E ^2u{J^) 
.1=0,1 



(4) 



J=0,1 



The corresponding total (3^^ transition strengths are: 



54 



E 



EA±(pn;aJ)Wj 



OJJ 
pn 



pn 



(5) 



For the discussion of the numerical results the following unperturbed moments are also 
needed: 



mlin = E (-)"A^(pn)A° (pn) (iVp^/^)' 7^^(pnpn; 



(6) 



pnL5 
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and 



mt{r) = (-)-^EA^(pn)A° (pn) {wH') /u^j,, (7) 

pn 

where uJ^^j-w are the unperturbed pn-energies for a given J'^, and 

A° (pn) = sJPpPn UpVn , A° (pn) = ^/p^u^v^. 

Finally, when the RHC mechanism is not considered, the half lives for the 0^ 0"*" 
transitions read 

where ^ is a kinematical factor pSI, is the nuclear matrix element, and the values of 



T are 



1 for [3132^ , 
^ for /?/5o. . 



3 Numerical Results and Discussion 
3.1 Unperturbed Matrix Elements 



As in our previous studies of the /3/3-decays [p, [10], |Tl], 0, the numerical calculations 
are performed with a 5-force (in units of MeV fm?) 

with different strength constants Vs and Vt for the particle-hole, particle-particle and pairing 
channels. 

For the nuclei ''^Ge, ^^5*6, ^^^Te, ^^^Te and ^°°Mo, we work in an eleven dimensional 
model space including all the single particle orbitals of oscillator shells ?)huj and ATiuj plus the 
0/i9/2 and O/iii/2 orbitals from the ^Tiu oscillator shell. The single particle energies (s.p.e.), 
as well as the parameters ff""'(p) and w|'"*'^(n), have been fixed by the procedure employed 
in ref. |]10| (i.e., by fitting the experimental pairing gaps to a Wood-Saxon potential well). 



8 



In fact, for the nuclei ^^Ge, ^"^Se, ^'^^Te, and ^^°Te we use here the same parameterization 



as in ref. 10 



For the nucleus ^^Ca the calculations have been done with three different model spaces, 
namely those including all the orbitals in the major shells: 2huj and 3huj {space A), Ohu to 
3hu) [space B), and Ohuj to 4huj [space C). Experimental single particle energies have been 
used for the orbitals, lpi/2, O/5/2, 1^3/2, O/7/2, lsi/2 and 0^3/2, here. For the remaining 
orbitals a single-particle energy spacing of huj = AlA^^^^ MeV is assumed. 

The unperturbed 2v and moments, given by eqs. @ and (|^) and displayed in table 
0, represent the upper limit to the perturbed moments. We discuss these first, in order to 
gain an insight into both the magnitude of (3/3 moments and the role played by the size of 
the configuration space. The Ou matrix elements have been calculated by taking 5 MeV as 
the mean excitation energy Uaj^ in the neutrino potential (|l|). It has been found that the 
largest fraction of the Oz/ strengths for all the nuclei is concentrated around this energy. 

We may note that, while the m2y[J'^) matrix elements are practically unchanged as 
we go to larger configuration spaces, the mQ^[J'^) matrix elements, and in particular those 
with J > 2, vary significantly. The total moment Mq^, on the other hand, shows only a 
small variation, due to the dominance of the mQ^[0^) and mQ^(l"'") in comparison with the 
moments mQ^[J > 2). This result is not surprising as, from the study of the charge-exchange 
resonances, we already know that for ^^Ca the spaces A and B are complete spaces only for 
the Fermi and GT transitions. For the first-forbidden resonances [L = 1,J'^ = 0^, 1^,2^), 
a calculation that does not include the Ahu shell as well, would be totally unacceptable. 
In the same way for the evaluation of the second-forbidden resonances with L = 2 and 
J'^ = l+,2+,3^ one should include the orbitals of the 5huj oscillator shell (otherwise the 
corresponding charge exchange sum rule will not be satisfied), and so on. Fortunately, the 
problem with the size of the configuration space is not so serious for the /?/5-decay. The 
reason is that the transition amplitude for this process is proportional to the factor UpV^UnVp 
that rapidly decreases as we move away from the valence orbitals. The same argument is not 
valid for charge-exchange transitions whose amplitudes carry the factor UpV^ or the factor 

UnVp- 
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3.2 Perturbed Matrix Elements 

The particle-hole channel parameter values, uf^ = 27 and vf^ = 64 for ^^Ca and uf'^ = 55 and 
^ph _ g2 ^YiQ remaining six nuclei, have been taken from a study of energy systematics 
of the GT resonances (|^9|, |6^. For further discussion it is convenient to introduce the 



parameters s and t, defined as the ratios between the T = 1, 5* = and T = 0, S = 1 
coupling constants in the PP channels and the pair ing force constants, i.e., 

2vPP 2vr 
" - t 



For a value of s = 1 the isospin symmetry is restored within the QRPA, leading to a 
concentration of 5+(0'^) strength in a single state and resulting in 

5_(0+) = ; m2^(0+) = 0, and mo^(0+) = 0. (8) 

From fig. |l] it is evident that the conditions (P) are fulfilled reasonably well for all the 
nuclei discussed here. As such we fix the spin-singlet PP strength at the value s = 1, 
obtaining m2u{0~^) = and A42u — i^2u{^^)- It should be noted that in most of the QRPA 
calculations performed so far, the moment m2u{0~^) has been simply ignored by invoking the 
isospin symmetry. Simultaneously, however, sizable values for the moment mou{0~^) have 
been reported. Apparently such calculations are inconsistent since the restoration of the 
isospin symmetry within the QRPA takes place at the level of the residual interaction and 
not at the level of the one body operator t+. 

Before looking at the experimental data we consider the question of validity of the BCS 
approximation for far off orbitals as the configuration space is extended. With this in mind, 
three different calculations have been performed for ^^Ca, namely 

Calculation I: Both the BCS and the QRPA equations have been solved within the single- 
particle space A. 

Calculation 11: The configuration space for solving QRPA equations extends over the major 
shells Qhjj — 4fiu (space C), but the gap equations have been solved within the space B only. 
The Ahuo shell is assumed to be totally empty. 

Calculation III: The full space C has been used at all steps of the calculation. 
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In each case, the pairing interaction strengths have been determined by fitting the experi- 
mental odd-even mass differences for the relevant nuclei. The calculated iS_(l"^), —m2^{l~^), 
— mou{^~^) and —Aio^, as a function of parameter t, are shown in fig. |^. The curves have 
been drawn up to the first pole of m2jy(l^), where the energy of the lowest virtual l"*" state 
becomes equal to the energy of the initial or final state and the QRPA breaks down. The 
magnitudes of all the quantities mentioned above depend to some extent on the size of the 
configuration space and the way in which the BCS equations have been solved. Yet, the gen- 
eral trend is the same in the three spaces. As a matter of fact, if one fixes the value of the 
parameter t at the lowest value of iS_(l"'"), we get practically the same results from all three 
calculations; that is: m2^(l+) = 0.09 MeV-\ mo^(l+) = 0.9 MeV'^ and Mo^ = -0.45 
MeV-\ 

We may generalize the conclusions drawn for the case of ^^Ca. That is, an enlargement 
of the configuration space, beyond two major oscillator shells, is not expected to modify in 
essence the model predictions. As such the 11 single-particle space described earlier should 
be a rather complete space for the description of the f3f3 processes in the remaining five 
nuclei. Fig. ^ shows the relevant results for the nuclei ^Ge, ^^5*6, ^^^Mo, ^^^Te and ^^°Te. 
From an inspection of figs. ^ and |^ the following common points may be made: 

1) The minimum of iS_(l'*") always occurs at a value of t = tsym that lies between the zero 
and the pole of m2u{i^)- 

2) As pointed out earlier in a similar context ([0, ^ ||, |, ||, |^), in the vicinity of tgym, the 
moment nioui^'^) is very sensitive to small variations of the parameter t, while the total 
moment Aiou is only moderately influenced by the same variation. The only exception is 
the nucleus '^^Ca. Here, due to a very pronounced dominance of the r?7oj,(l"'") moment over 
the moments with higher multipolarities, the total Ou matrix element passes through zero 
at a value of t that is very close to tgym- 

3.3 Comparison with Experimental Data 

The measured 2u and Ou half-lives are listed in table |] along with the values of kinemat- 
ical factors G'^" and G^'^, appropriately renormalized for an effective axial- vector coupling 
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constant qa = —gv, and the resulting observables \J^2u\ and lA^Oi/ ^™"^ I- 

The results of the calculations for the matrix elements A^2i/ and Aiou, as well as for 
the predicted 2z/ half-lives and the neutrino masses, are shown in table ^ From the upper 
panel of this table it is seen that when t = tgym the calculated 2v matrix element for '^^Ca 
does not contradict the experimental limit and that the measurement of the moment M.2u 
in ^'^Se is well accounted for by the theory. But, the calculated matrix elements J^2u turn 
out to be too small for ^Ge and ^^^Mo and too large for ^^^Te and ^^^Te (in both the cases 
by a factor of ~ 3). It is worth noting that the measured values of |A^2!/| for "^^Ge and 
^^^Mo nuclei are not very different from the calculated ones even when no PP interaction 
is included (see fig. From the same figure, we may notice that the minima of iS_(l"'") 
are not so well defined and it is precisely near these minima that the calculated values of 
M.2u vary rather abruptly. Evidently, as such, the experimental data for the 2v half-lives 
can be reproduced in all the six nuclei with a value of t very close io t = tsym- Finally, the 
calculated M.qv moments (and estimated neutrino mass limits) for t- values that reproduce 
the measured matrix elements M.2v, when these are assumed to be positive (t = t^) and 
negative (t = t^), are listed in middle and lower panels of table ^, respectively. One notices 
that in all the cases tgym — t^, and tgym and lead to practically the same upper limits for 
effective neutrino mass <my>= leV. 



4 Conclusions 

We have calculated the 2i/ and Oz/ double beta observables for the nuclei ^^Ca, "^^Ge, ^"^Se, 
^™Mo, ^^^Te and ^^°Te in the framework of the QRPA model satisfying the constraints 
imposed on the particle-particle coupling strengths by the isospin and SU(4) symmetries, 
i.e., s = 1 and t = t^ym- With the parameter t equal to tgym the calculations are consistent 
with the experimental limit for the 2z/ lifetime of ^^Ca and they correctly reproduce the 2i> 
lifetime of ^'^Se. For the remaining nuclei the 2i/ half lives are either overestimated (for "^^Ge 
and ^°°Mo) or underestimated (for ^^^Te and ^^°Te) by an order of magnitude. This does not 
imply that one has to forgo the idea of reconstructing the SU(4) symmetry by the residual 
interaction. We believe that the restoration of both the isospin and SU(4) symmetries is a 
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genuinely useful feature of QRPA and undoubtedly plays an important role in the intricate 
physics involved in jS/S processes. One also should bear in mind that: i) with a comparatively 
small variation ( < 10%) of t with respect to tsym, i-e., with t = t-^ it is possible to account 
for the in the cases, and ii) the minimum value of the GT 5_(1+) strength critically 
depends on the spin-orbit splitting over which we still do not have a complete control. In the 
same context it should be interesting to study the role of forbidden virtual states, particularly 
those with J'^ = 1~ and 2~ (see table ^), in building up the total J^2u moments . 

For a long time we have been worrying about the completeness of the virtual states in the 
evaluation of the Ou moments, or in other words: how the Oz/ moments depend on the size of 
the configuration space? The results displayed in table || and fig. || seem to suggest that the 
enlargement of the space beyond two oscillator shells does not have a significant effect either 
on the GT strength iS_(l"'') or on the f3(3 moments 'm2iy{l~^), mo;y(l"^) and A4ou- In the first 
case the configuration space turns out to be sufficiently complete while in the second case 
the pairing factor dependence of the (3f3 transition amplitudes inhibits the far off orbitals to 
contribute. We feel, however, that this question is not yet resolved unambiguously. 

A final and rather general comment is in order. Besides the issue of the procedure adopted 
for fixing the particle-particle strength parameter, there are some additional problems within 
the QRPA calculations of the matrix element Ai2u, as yet not fully understood. They are 
related with the type of force, choice of the single particle spectra, treatment of the difference 
between the initial and final nuclei, etc. All these things are to some extent uncertain and 
therefore it is open to question whether it is possible, at the present time, to obtain a more 
reliable theoretical estimate for the 2i/ half lives that the one reported here. Similar remarks 
stand for the A4ou moments and hence for the neutrino mass limits. The difference in a 
factor of about 2 — 3 between both: i) the results obtained by the Pasadena group and 
the groups of Tiibingen and Heidelberg for "^^Ge and ^"^Se nuclei, and ii) the previous and 
present calculations for ^^^Mo, ^^^Te and ^^°Te nuclei, is just a reflection of the unavoidable 
uncertainty of the QRPA calculations, and it is difficult to assess which one is "better" and 
which is "worse". 
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A Radial Matrix Elements for the Oz//3/3-decay 



Within the Horie-Sasaki formahsm the radial matrix elements (|]) for the harmonic 
oscillator wave functions read as 

7^^(pnpV;o;„J.) = R [M(p, n)M(p', n')]^^ 



X 



am(p, n)a„/(p', n')f^{m, m'; Uaj^), 



with 



M(ni, n'i') = 2"+"'n!n'!(2^ + 2n + l)!!(2f + 2n + 1)!!, 

i'\ (2£ + 2n+ 1)!! (2f + 2n' + 1)!! 

\ k l \ I 

{k+k'=s) 



(k+k'=s) / 



{2i + 2k + l)\\ {2i' + 2k' + 



„j , . , X ^ (m~L m' — L \ 



and 



where v = Mu/h is the oscillator parameter. For v(g; uJaj^) given by eq. (|T]), Tomoda et al. 



5^ have obtained the following recurrence relation for the momentum space integrals 

2^(2/i-l)!! V2^ 



TT 2/^ 



TT 



-u(/i- l)! + uV^-i(u) 



Jo(u) = W u$(u) 

V TT TT 



where u = uJaj^ / \/2ij and 



$(u) 



(it = ea;p(— u^) 



roo \ 

v^TT / exp{t^)dt - -Eiiy?) 
Jo 2 



10 t + n 

When finite nucleon size (FNS) effect and the short range (SR) two-nucleon correlations 
are included, the potential v{q;Uaj^) takes the form 



^ENs+sRiq; ^aj^) = ^FNsiq; ^aj^) - Av(g) + A'v(g), 
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with 



Av'(g) 



TT 



+ q^J 



, Av(g) = —In 



q + qc 



.n=l 



Q-Qc 
A2 



A.^ + iq±qc 



where A = 850MeV is the cutoff for the dipole form factor and q = 3.93 fm^^ is roughly 
the Compton wavelength of the a;-meson. The corresponding integrals ^^{uJaj^) have to be 
evaluated numerically. 
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Figure Captions 



Figure 1: Fermi observables S^(J'^ — 0+), mzuiJ^ — 0^) (in units of [MeV^] ^) and mo^{J'^ — 
0+) for the nuclei ^«Ca, ^^Ge, ^^^e, i°°Mo, ^^^Te and ^^^Te, as a function of particle-particle 
S — 0, T — 1 coupling constant s. 

Figure 2: Gamow- Teller observables S-{J'^ = l"*"), m2v{J^ = l"*") (in units of [Me\^]~^), 
mQu{J'^ = 1+) and the total Aiou moment for the nucleus '^^Ca, as a function of the particle- 
particle S = 1, T = coupling constant t, within the single-particle spaces A, B and C. 

Figure 3: Gamow- Teller observables S-{J'^ = 1"*"), mzv^J^ = l"*") (in units of [Mel^]~^), 
mo^(J" = 1+) and the total Mou moment for the nuclei ^^Ge, ^^Se, ^^^Mo, ^^^Te and ^^°Te, 
as a function of the particle-particle S — 1, T — coupling constant t. 
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Tables 



Table 1: Upper bounds on the effective neutrino mass < rriu > (in eV) obtained from the 
QRPA calculations of the nuclear matrix elements. For the sake of comparison, in all the 
cases the effective axial vector coupling constant g"^^ = —gy has been employed. This means 
that the results for < m,^ > from refs. 0, [T^, ^ have been properly renormalized by the 
factor (1.25)^. The experimental data for the half-lives which have been used are indicated 
in the second row. The two values of the Pasadena group correspond to their PP strengths: 
(a) a[ = -390 and (b) a[ = -432, in units of MeV fm^. 



Exp. 



ref. 
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ref. m 



ref. m 



ref. 



2|] 



ref. m 



ref. [gS 



ref. 
ref. 
ref. 
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(aj 
(b) 



2.3 
4.4 
10 
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20 
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20 
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2.4 
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1.5 



24 
22 
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21 
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m 




3.1 


12 


3.8 


31 
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Table 2: Unperturbed m2u{J^) and niQuiJ^) moments in units of [Mel^] ^. As explained in 
the text, three different single-particle spaces have been used for the nucleus ^®Ca. 























A 


B 


C 












""Ivy" J 

0+ 


0.249 


0.255 


0.253 


0.423 


0.491 


0.256 


0.678 


0.615 


1+ 


0.406 


0.417 


0.414 


0.923 


0.966 


1.513 


1.602 


1.434 


0+ 


0.928 


0.953 


0.990 


2.258 


2.424 


1.857 


2.671 


2.472 


1+ 


2.072 


2.158 


2.268 


6.168 


6.259 


7.204 


7.310 


6.751 


2+ 


0.316 


0.384 


0.426 


1.285 


1.283 


1.470 


1.483 


1.367 


3+ 


0.321 


0.394 


0.438 


1.112 


1.179 


1.421 


1.260 


1.164 


4+ 


0.109 


0.131 


0.151 


0.507 


0.525 


0.708 


0.636 


0.587 


5+ 


0.134 


0.154 


0.170 


0.407 


0.459 


0.636 


0.497 


0.464 


6+ 


0.037 


0.037 


0.044 


0.196 


0.214 


0.356 


0.253 


0.236 


7+ 


0.075 


0.075 


0.079 


0.145 


0.187 


0.275 


0.212 


0.200 


8+ 


0.000 


0.000 


0.001 


0.055 


0.070 


0.160 


0.096 


0.092 


9+ 


0.000 


0.000 


0.000 


0.056 


0.083 


0.065 


0.087 


0.084 


10+ 


0.000 


0.000 


0.000 


0.002 


0.002 


0.006 


0.029 


0.028 


0- 


0.027 


0.040 


0.071 


0.177 


0.194 


0.337 


0.130 


0.126 


1- 


0.307 


0.378 


0.564 


1.604 


1.749 


2.655 


1.367 


1.315 


2- 


0.209 


0.250 


0.337 


1.161 


1.266 


1.628 


1.128 


1.078 


3- 


0.158 


0.183 


0.244 


0.882 


0.954 


1.124 


0.929 


0.881 


4- 


0.061 


0.061 


0.093 


0.517 


0.559 


0.735 


0.559 


0.531 


5- 


0.062 


0.063 


0.087 


0.417 


0.453 


0.456 


0.479 


0.453 


6- 


0.000 


0.000 


0.010 


0.157 


0.175 


0.276 


0.249 


0.237 


7- 


0.000 


0.000 


0.009 


0.181 


0.197 


0.171 


0.258 


0.245 


8- 


0.000 


0.000 


0.000 


0.014 


0.019 


0.051 


0.084 


0.082 


9- 


0.000 


0.000 


0.000 


0.015 


0.018 


0.048 


0.128 


0.122 


-Ml 


4.816 


5.261 


5.983 


17.317 


18.268 


21.641 


19.844 


18.516 
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Table 3: Values of the measured half-lives T1/2, kinematical factors G and the experimental 
A^F-quantities for the 2z/ and Ou decays. The G factors are those from ref. but 
renormalized for qa = —Qv- 



Observable ^^Ca ^ ^Gt ^'ge ^°°Mo ^^^Te ^^""Tt 

1+u.y b 1 nc+u.^ti c n 1 1 c+u.uau d '7 _i_ n a\ in4e 



Tyjyr 10^0] >0.36'^ 9.2^^:1'' l-OST^ioe 0.1151^^^ (7-7 ± 0.4) 10^ 27 ± 1 " 

G2'^[?/r(Mel^)2]-i 0.423 lO^^^ 0.139 lO'^^ 0.464 lO'^^ 0.101 10"^^ 0.911 lO'^^ 0.512 10"^^ 

\M^,\\MeVY^ < 0.081 0.280t°:°?^ 0.14ll°:0?l 0.294l°:01 0.0381°;°} 0.0271°;°} 

Tll^\yrW^\ >9.5-^ > 1400^^ > 27'= > 4.7"^ > 7700^ >2.5'^ 

G°'^[l/r]-i 0.260 10-13 0.261 10"^^ 0.114 lO'^^ 0.187 10-^3 0.746 lO^^^ 0.18110-^3 

lA^o^^^l 10^ <0.64 <0.17 <0.57 < 1.07 < 0.13 < 1.5 



(laboratory data) ref. [^| 



(laboratory data) ref. [32| 
(laboratory data) ref. p5| 
(laboratory data) ref. [p^ 
(geochemical data) ref. |^ 
(laboratory data) ref. |2^ 
(laboratory data) ref. 
(laboratory data) ref. 
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Table 4: Calculated 2v and Oz/ moments and the corresponding upper limits for the effective 
neutrino mass < my >. Here s = f and three different sets of parameter t have been 
considered, namely, t = tsym, t = and t = ti- The last two reproduce the measured 2i/ 
matrix element (without taking error bars into consideration) when these are assumed to be 
positive and negative, respectively. For t — tsym, the calculated T^/'2 values are also shown. 
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-7.0 
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71 
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1.0 


12 




0.28 


71 
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0.95 


7.9 X 10=^ 


2.1 




1.48 
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1.39 


1.53 


1.35 


1.30 




-0.48 


-4.5 


-5.3 


-6.9 


-7.3 


-7.0 


{m^)^ [eV] 


68 


1.9 


7.1 


10 


1.0 


11 


h 


1.20 


0.30 


1.00 


1.30 


1.22 


1.20 




-1.41 


-10.2 


-7.7 


-8.3 


-8.1 


-7.5 


(m.) , [eV] 


23 


0.87 


4.2 


6.6 


0.86 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9312007vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9312007vl 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9312007vl 



